ON THE TWO-DIMENSIONAL MOMENT PROBLEM 
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Abstract. In this paper we obtain an algorithm towards solving the two- 
dimensional moment problem. This algorithm gives the necessary and sufficient 
conditions for the solvability of the moment problem. It is shown that all 
solutions of the moment problem can be constructed using this algorithm. In a 
consequence, analogous results are obtained for the complex moment problem. 



1. Introduction and preliminaries 

In this paper we analyze the two-dimensional moment problem. Recall that 
this problem consists of finding a non-negative Borel measure /i in R 2 such that 

x^x^dfi = s mt n, m,n G Z + , (1.1) 

where {s m , n }m,nez + is a prescribed sequence of complex numbers. 

The two-dimensional moment problem and the (closely related to this subject) 
complex moment problem have an extensive literature, see books [7], [1], [3], sur- 
veys [5], [4] and [8]. Some conditions of solvability for this moment problem were 
obtained by Kilpi and by Stochel and Szafraniec, see e.g. [1] and [8]. However, 
these conditions are hard to check. Putinar and Vasilescu derived conditions of 
solvability and a description of all solutions by means of a dimensional exten- 
sion [6] (even for the iV-dimensional moment problem). The two-dimensional 
moment problem is solvable if and only if the prescribed sequence of moments 
can be extended to a sequence {s mt n t k}m,n,kez + , satisfying some easy conditions 
(including the positivity condition). This extended sequence is the moment se- 
quence for an extended moment problem: 

/ ^T x 2~, 7 rrrdyu = s min>fc , m,n,keZ + . (1.2) 

Jvt? (1 + x\ + xl) 

The unique solution of the moment problem (1.2) provides a solution of the 
two-dimensional moment problem. In this way all different extensions define all 
different solutions of the two-dimensional moment problem. However, it is not 
clear whether such extensions exist and what is a procedure for the construction 
of extensions {s m ^ )k } mtnik&+ - 

The method of our investigation uses an abstract operator approach, see [9]. 
Firstly, we obtain a solvability criterion for an auxiliary extended two-dimensional 
moment problem. This moment problem is somewhat similar to the moment 
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problem (1.2) but we do not see any direct relationship. It is shown that the 
extended two-dimensional moment problem is always determinate and its solution 
can be constructed explicitly. 

An idea of our algorithm is to extend the symmetric operators related to the 
two-dimensional moment problem, not "entirely", but on a discrete set of points. 
It is shown that all solutions of the moment problem (1.1) can be constructed on 
this way. Roughly speaking, the final algorithm reduces to the solving of finite 
and infinite linear systems of equations with parameters. 

In a consequence, analogous results are obtained for the complex moment prob- 
lem. 

Notations. As usual, we denote by M, C, N, Z, Z + the sets of real numbers, 
complex numbers, positive integers, integers and non-negative integers, respec- 
tively. The real plane will be denoted by M. 2 . For a subset S of M 2 we denote by 
%$(S) the set of all Borel subsets of S. Everywhere in this paper, all Hilbert spaces 
are assumed to be separable. By (•, •)# and \\-\\h we denote the scalar product and 
the norm in a Hilbert space H, respectively. The indices may be omitted in obvi- 
ous cases. For a set M in H, by M we mean the closure of M in the norm || • \\h- 
For {a; m ,„} m ,„ e z + , x m ^ n G H, we write Lin{x mjn } mjneZ+ for the set of linear combi- 
nations of elements {x m>n } mineZ+ and span{x m>n } m>neZ+ = Lm{x m>n } m , n£Z+ . The 
identity operator in H is denoted by Eh- For an arbitrary linear operator A in 
H, the operators A*, A, A -1 mean its adjoint operator, its closure and its inverse 
(if they exist). By D(A) and R(A) we mean the domain and the range of the op- 
erator A. The norm of a bounded operator A is denoted by \\A\\. By Pff = Ph 1 
we mean the operator of orthogonal projection in H on a subspace Hi in H . By 
L 2 we denote the usual space of square-integrable complex functions f(xi,X2), 
x±,X2 G M 2 , with respect to the Borel measure U, in M 2 . 

2. The solution of an extended two-dimensional moment problem. 

Consider the following moment problem: to find a non-negative Borel measure 
[i in R 2 such that 

x™(xi + i) k {xi - %fxV 1 {x2 + i) r {x 2 - ifdfi = u mik ,i; n , r ,t, 

m,n<E'Z + , k,l,r,t&Z, (2.1) 

where {u mj k,i;n,r,t}m,nez+,k,i,r,tez is a prescribed sequence of complex numbers. 
This problem is said to be the extended two-dimensional moment problem. 

We set 

Q = {(to, k, l; n, r,t) : m,n G Z+, k, l,r,t G Z}, 
Q = {(m, k, l; n, r,t) : m,n G Z + , k, I, r, t G Z, k = I = r = t = 0}, 

n' = n\n . 

Let the moment problem (2.1) have a solution \i. Choose an arbitrary function 

P(x u x 2 ) = ^ a m,k,i;n,r,txT{xi + i) k (xi - i) l x%(x 2 + i) r (x 2 - i)*, 

(m,k,l;n,r,t)£Q 
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where all but finite number of complex coefficients a m ,k,i;n,r,t are zeros. Then 



0< / \P{x u x 2 )\ 2 dfJt = Yl 



^m,k,l;n,r,t0^m' ,k' ,l';n' ,r' ,t' 

(m,k,l;n,r,t),(m' ,k' ,l';n' ,r' ,t')en 

< +m '(xi + i) k+l '(x 1 - i) l+k ' x% +n ' (x 2 + i) r+t '(x 2 - i) t+r 'dix 



— ^ a m,k,l;n,r,t a m' ,k' ,l';n' ,r' ,t' u m+m' ,k+V \l+k';n+ri ,r+t' ,t+r' ■ 

(m,k,l;n,r,t),(m' ,k' ,V ;n' ,r' ,t')eU 

Therefore 

^ ^ Oim,k,l\n,r,t a m' ,k' ,V\n' ,r' ,t' u m+m' ,k+l' ,l+k';n+ri ,r+t' ,t+r' ^ 0, 

(m,k,l;n,r,t),(m',k',l';n',r',t')&l 

(2.2) 

for arbitrary complex coefficients a m ,k,i;n,r,ti where all but finite number of a m ,k,i;n,r,t 
are zeros. The latter condition on the coefficients a m ,k,i,n,r,t in infinite sums will 
be assumed in similar situations. 

We shall use the following important fact (e.g. [2, pp. 361-363]). 

Theorem 2.1. Let a sequence of complex numbers {u mt k,i- n ,r,t}(m,k,i;n,r,t)en satisfy 
condition (2.2). Then there exist a separable Hilbert space H with a scalar product 
(■,-)h and a sequence {x m ^i. n ^ t }( m ,k,i;n,r,t)^n in H, such that 

(Xm,k,l;n,r,ti x m' ,k' ,V \n' ,r' ,t') H = U m ^- m i t k+l' ,l+k';n+n' ,r+t' ,t+r' , 

(m, k, I; n, r, t), {m! , k', I'; n\ r', t') e fi, (2.3) 

and Spa.n{x mj k,l-n,r,t}(m,k,l;n,r,t)Gn — H. 

Proof. (We do not claim the originality of the idea of this proof). Choose an arbi- 
trary infinite-dimensional linear vector space V (for instance, one may choose the 
space of all complex sequences (u n ) ne ®, u n e C). Let X = {x m ,k,i;n,r,t}(m,k,i;n,r,t)en 
be an arbitrary infinite sequence of linear independent elements in V which is in- 
dexed by elements of Q. Set L x = ^{x mik ,i ; n,r,t}{m,k,i;n,r,t)en- Introduce the 
following functional: 

l x ,y]= Ctm,k,l;n,r,t{3m',k',l';ri,r',t' 
(m,k,l;n,r,t),(m' ,k' ,l';ri ,r' ,t')en 

* u m+m',k+l',l+k';n+n',r+t',t+r' > (2-4) 

for x,y e L x , 

x ^ CV-in,k,i,n,r,t x m,k,i,n,r,t] 

(m,k,l;n,r,t)£tt 

y = @m> \k> \l <;n' >' ,t'Xm> ,k' ' ,l> ';n' >' ',t> 
(m',k',l';n',r',t')en 

where a m , fc ,/ ; n,r,t, Pm',k',i'-,n>,r',t' G C. Here all but finite number of indices a m , fc ,z ; n,r,t, /3 m >,k>,i>;n',r>,t> 
are zeros. 

The set Lx with [•, •] will be a quasi-Hilbert space. Factorizing and making the 
completion we obtain the desired space H (e.g. [3]). □ 
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Let the moment problem (2.1) be given and the condition (2.2) hold. By The- 
orem 2.1 there exist a Hilbert space H and a sequence {x m ^,i- n ,r,t}{m,k,i;n,r,t)^j * n 
H, such that relation (2.3) holds. Set L = Lm{x mt k,i- n ,r,t}(m,k,i,n,r,t)en- Introduce 
the following operators 

-Aq ^ ^ &m,k,l;n,r,t%m,k,i,n,r,t ^ ^ ^m,k,l;n,r,tXm+l,k,l;n,r,ty (^'5) 

(m,k,l;n,r,t)£tt (m,fc,/;n,r,t)Sf2 

Bq ^ Qm,k,l;n,r,t%m,k,l;n,r,t ^ Qm,k,l;n,r,t%m,k,l;n+l,r,ti (^"^) 

(m,k,l;n,r,t)sQ (m,k,l;n,r,t)GQ 

where all but finite number of complex coefficients a mi k,i;n,r,t are zeros. Let us 
check that these definitions are correct. Indeed, suppose that 

(m,k.l;n,r,t)en (m',k',l';n',r',t')€U 

(2.7) 

We may write 



^ ^ ^ / -m,k,l;n,T,tXm+l,k,l;n,r,t: X a ,b,c;d,e,f 

(m,k,l;n,r,t)££l 



H 



(m,i;,!;ii,r,t)en 

= £ ^W;»,r,*(^ A «,r,*,^ 1AeiW )H 

(m,,fc,i;ra,r,i)Gn 

= (x, x a+ i ;6iC;diei /)ij-, (a, b, c; d, e, f) G fi. 
In the same manner we obtain: 



a m',k',l';n',r',t' x rn'+l,k',l';n',r',t' , x a,b,c;d,e,f 

(m',k',l';n',r',t')en 
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= (x, x a+ltb>c . dte j) H , (a, b, c; d, e, f) G fl 
Therefore the definition of Aq is correct. The correctness of the definition of Bq 
can be checked in a similar manner. Notice that 

{■^■0Xm,k,l;n,r,t j •£a,6,e;d,e,/)iT (^m+l,fc,l;n,r|f > x a,b,c;d,e,f)ll 
^m+l+a,fc+c,£+&;?i+c(,r+/,t+e {Xm,k,l\n,r,t ; Xa+l,6,c;d ; e,/)iy 

= (^,t,i in ,r,t, A)X a ,b, c; d,e,/)i?, (m, fc, i; 7i, r, (a, 6, c; d, e, /) G fi. 
Therefore A is symmetric. The same argument implies that Bq is symmetric, as 
well. 

Suppose that the following conditions hold: 

u m+l+a,k+c,l+b;n+d,r+f,t+e + lu m+a,k+c,l+b;n+d,r+f,t+e 

= u m+a.k+l+c,l+b;n+d.r+f,t+e, (2-8) 
u m+l+a,k+c,l+b;n+d,r+f,t+e ~ tu m+a,k+c,l+b;n+d,r+f,t+e 
r U j m+a,k+c,l+l J rb\n J rd,r+ f ,t+ei 

(2.9) 
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u m+a,k+c,l+b;n+l+d,r+f,t+e + iu m+a,k+c,l+b;n+d,r+f,t+e 
= u m+a,k+c,l+b;n+d,r+l+f,t+ey 

(2.10) 

u m+a,k+c,l+b;n+l+d,r+f,t+e ~ lu m+a,k+c,l+b;n+d,r+f,t+e 

= u m+a,k+c,l+b;n+d,r+f,t+l+e> (2-11) 

for all (m, k, I; n,r, t), (a, b, c; d, e, /) G fl These conditions are equivalent to 
conditions 

(■Kin+l,k,l;n,r,t '^•^m,k,l;n,r,ti •^a,b,c;d,e,,f) H 

= (Xm,k+l,l;n,r,t,Xa,b,c;d,e,f)H, (2-12) 
\%m-\-l,k,l;n,r,t ^■£m,k,l;n,r,ti 2"a,fe,c;d,e,/) H 

\%m,k,l+l;n,r,ti x a,b,C]d,e, f) H , (2.13) 
\3>m,k,l;n+l,r,t ^•^m,k ! l;n,r,tj •^a,b,c;a!,e,/) H 

\%m,k,l;n,r+l,ti x a,b,c;d,e, f)n, (2.14) 
(■^m,k,l;n-\-l,r,t i x m,k,l;ri,r,ti ■%a,b,c;d,e,f*) H 

\%rn,k,l;n,r,t+l) x a,b,c;d,e, fh, (2.15) 

for all (m, k, I; n, r, t), (a, b, c; e, /) G f2. The latter conditions are equivalent to 
the following conditions: 

^m+l,k,l;n,r,t ^m,k,l;n,r,t •^"m,k+X,l;n,r,ti (2-16) 

■^m+l,k,l;n,r,t ^m,k,l;n,r,t ■^m,k,l+l;n,r,t i (2-17) 

■Em,k,l;n+l,r,t ^■^m,k,l;n,r,t •^m,k,l;n,r+l,t ; (2-18) 

•Em,k,l;n+l,r,t ^m,k,l;n,r,t %m,k,l;n,r,t+l ; (2-19) 

for all (m, k, I; n, r, t) G f2. The last conditions mean that 

(Aq + iE H )x m ^J- n ,r,t = ^m,fc+l,i;n,r,t; (2.20) 

(A — iE}j)x m ^k.l-n,r,t — x m,k,l+l;n,r,t, (2-21) 

(Bq + iE H )x m ^,l- n ,r,t = Xrn,k,l;n,r+l,ti (2.22) 

(Bq — iEH)x mt k,l ]n ,r,t — Xm,k,l;n,r,t+1, (2.23) 

for all (m, k, I; n, r, t) G Q. The latter conditions imply that 

(A ±iE H )L = L, (B ±iE H )L = L. 

Therefore operators Aq and Bq are essentially self-adjoint. The conditions (2.20)- 
(2.23) also imply that 

(Aq + iE H ) 1 X m ^,l;n,r,t = ^m,fc-l,«;n,r,t , (2.24) 

(Aq — lEu) l X m ,k,l;n,r,t = Xm,k,l-l;n,r,ti (2.25) 

(Bq + iE H ) 1 X mt k t l ;n ,r,t — Xm,k,l;n,r-l,t, (2.26) 

(Bq — iEff) 1 X m ,k,l;n,r,t = ^m,fe,«;n,r,t-l, (2.27) 

for all (m, k, /; n, r, t) G Q. 

Consider the Cayley transformations of Aq and Bq: 

V Ao = (A - iE h )(Aq + iE H )~\ (2.28) 

V Bo = (Bq - iE h )(Bq + iE H y\ D(Aq) = D(Bq) = L. (2.29) 
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(2.30) 



By continuity we extend the isometric operators Va and V# to unitary operators 
Ua and Vb q in H, respectively. By continuity we conclude that 



Set A = Aq, B = Bq. The Cayley transformations of the self-adjoint operrators 
A and B coincide on L with Ua and Ub , respectively. Thus, the Cayley trans- 
formations of A and B are Ua and Ub , respectively. Therefore, operators A and 
B commute. 
Notice that 



x m ,k,i;n,r,t = A m (A + i)\A - i) l B n (B + i) r (B - *)'x ,o, ;0,o,o, (2-32) 



and then by substitution of each relation into previous one we obtain rela- 
tion (2.32). 

For the commuting self-adjoint operators A and B there exists an orthogonal 
operator spectral measure E{x) on <B(R 2 ) such that 



U Ao U Bo x = U Bo Ua x, 



x g H. 



(2.31) 



for all (m, k, I; n, r, t) G Q. In fact, by induction we can check that 

X m ,k,l;n,r,t = (B — iEtff 'X m ,k,l;n,r,0, t G Z, 

for any fixed m, n G Z + , k,l,r G Z; 

Xm,k,l;n,r,0 = (B + iE H ) T Xm^nfifi, T G Z, 

for any fixed m,n G Z + , fc,ieZ; 

£m,fc,Z;n,0,0 = ^ ^m,fc,t;0,0,0) n ^ Z + , 

for any fixed m G Z + , A;, / G Z; 

^^1,^,250,0,0 — (^4 — ^H) Z ^m,fc,0;0,0,0) ^ £ Z, 

for any fixed m G Z + , /c G Z; 

^m,fe,0;0,0,0 = (^ + «-E'i?)^m,0,0;0,0,0) A; G Z, 

for any fixed m G Z + ; 

^m,0,0;0,0,0 = ^0,0,0:0,0,0 > m £ Z + , 




(2.33) 



Then 





i*(xi + i) fc (xi - i)'x2 (x 2 + «) r (x 2 - z)*d(#(x)xo,o,o;0,o,o, ^o,o,o ; o,o,o)h- 
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Hence, the Borel measure 

M — (E (x)xofi fi-fi fifi, £o,0,0;0,0,o),tf) (2.34) 

is a solution of the moment problem (2.1). 

Theorem 2.2. Let the extended two-dimensional moment problem (2.1) be given. 
The moment problem has a solution if and only if conditions (2.2) and (2.8)- 
(2.11) are satisfied. If these conditions are satisfied then the solution of the mo- 
ment problem is unique and can be constructed by (2.34)- 

Proof. The sufficiency of conditions (2.2) and (2.8)-(2.11) for the existence of a 
solution of the moment problem (2.1) was shown before the statement of the 
Theorem. The necessity of condition (2.2) was proved, as well. Let us check 
that conditions (2.8)-(2.11) are necessary for the solvability of the moment prob- 
lem (2.1). 

Let fj, be a solution of the moment problem (2.1). Consider the space L 2 ^ and the 
following subsets in L^: 

L/i = Un{x^(x 1 + i) k (x 1 -i) l x^(x 2 + i) r (x 2 -i) t }(^ m ,k,i-,n,r,t)en, = V (2.35) 
We denote 

y m ,k,i;n,r,t ■= xT(xx+i) k (x x -i) 1 x%(x 2 +i) r (x 2 -if , (m, k, I; n, r, t) G ft. (2.36) 
Notice that 

{ym,k,l;n,r,tiVm',k',l';n',r',t')Lf l — ^m+m' ,k+l' ,l+k';n+n' ,r+t' ,t+r' i (2.37) 

for all (m, k, l; n, r, t), (m', k', I'; n', r', t') G f2. Consider the operators of multipli- 
cation by the independent variable in L^: 

A fM f(xi,x 2 ) = x 1 f(x 1 ,x 2 ), B fl f(xi,x 2 ) = x 2 f(x!,x 2 ), f G L 2 . (2.38) 
Notice that 

(An + iE L 2))y mikil . nirtt = y m> k+i,i;n,r,t, (2.39) 

(Afj, — iEu2jy mj k,l- n ,r,t — y m ,k,l+l;n,r,t, (2.40) 
(Bfj, + iEi2)y mk ^ nr ^ t = y m ,k,l;n,r+l,t, (2-41) 
(Bfj, — iE L 2)y m)k j. n)r . jt = y m ,k,l;n,r,t+l, (2.42) 

for all (m, k, I; n, r, t) , (m', k', l';n',r', t') G VL. 

Since conditions (2.2) are satisfied, by Theorem 2.1 there exist a Hilbert space H 
and a sequence of elements {x m ,k,i;n,r,t\ (m,k,i;n,r,t)&ii i n H, such that relation (2.3) 
holds. Repeating arguments after the Proof of Theorem 2.1 we construct opera- 
tors Aq and Bq in H. Consider the following operator: 

Wo Ot m ^k,l-n,r,tym,k,l;n,r,t = Qm,k,l-,n,r,t%m,k s l\n,T,ti (2.43) 

(m,k,l;n,r,t)£Q (m,k,l;n,r,t)£Q 

where all but finite number of complex coefficients a m ^,i- n ,r,t are zeros. Let us 
check that this operator is defined correctly. In fact, suppose that 

(m,k,l;n,r,t)&l (m',k',l';n',r',t')eU 

(2.44) 
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where f3 m ' k' i'-ri r' t' £ C. We may write 







^ ^ {.&rn,k,l;n,r,t fim,k,l;n,r,t)ym,k,l;n,r,t 
(m,k,l;n,r,t)£tt 



hi 



^ ^ (c^m,fc,i;n,r,t Pm,k,l:n,r,t) 

(m,k,l;n,r,t),(m' ,k' ,l';ri ,r' ,t')eU 



*(am',fc',/';n',r',i' ~ Pm> ,k> ,l';n' ,r' ,t')(ym,k,l;n,r,t, Urn' ,k' ,V ;n' ,r' ,t') 

^ ^ {.&rn,k,l;n,r,t Pm,k,i,n,r,t) 

(m,k,i,n,r,t),(m' ,k' ,V;n' ,r' ,t')£fl 



: (oim',k',l';n',r',t' ~ Pm' ,k' ,1' ;n' ,r' ,t')(x m ,k,l;n,r,t, %m' ,k> ,l';ri ',r> \f) 



II 



H 



^ ^ {p^m,k,l;n,r,t Pm,k,l;n,r,t)-^m,k,l;n,r,t 
(m,k,l;n,r,t)sQ 

Thus, the operator Wq is defined correctly. If x G H and 

^ ^fm,k,l;n,r,tym,k,l;n,r,ti 
(m,k,l;n,r,t)£Cl 

where ■y m ,k,i;n,r,t G C, then 
(W x,W x) H = 

0^m,k,l;n,r,t'ym' ,k' ,V ;n' ,r' ,t' 

(m,k,l;n,r,t),(m' ,k> ,l';ri ,r> ,t')eO 

*{Xm,k,l;n,r,h x m',k',l';n',r',t')H 

= 0£ mt k t l ;n ,r,tlm',k',l';n',r',t'(ym,k,l-,n,r,t, Vm' ,k' ,V \n> ,r> ,t') H 

{m,k,l;n,r,t),(m' ,k' ,V;n' ,r' ,t')&l 

= {x,x)i%. 

By continuity we extend Wo to a unitary operator W which maps if M onto H. 
Observe that 

W 1 AoWy m! k,l;n,r,t = y m +l,k,l;n,r,t = ^nl/m,fcJ;ri,r,(, (2-45) 
1 B Wy mj k ! l-n,r,t = y m ,k,l]n+l,r,t = -B/x?/m,fc,Z;n,r,t , (2.46) 

for all (m,k,l;n,r,t) G f2. By using the last relations in relations (2.39)-(2.42) 
we obtain relations (2.20)-(2.23). The latter relations are equivalent to condi- 
tions (2.8)-(2.11). 

Let us check that the solution of the moment problem is unique. Consider the 
following transformation 

T : (xi,x 2 ) G M 2 \-> (ip, VO G [0, 2tt) x [0, 2tt), 



up 



Xi + i 



x 2 + i 



and set 



X\ — i x 2 — i 

u(A) = //(T-^A)), A G <B([0,2tt) x [0,2tt)). 



(2.47) 
(2.48) 
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Since T is a bijective continuous transformation, then v is a non-negative Borel 
measure on [0, 2n) x [0,27r). Moreover, we have 



U0,k,-k:0,L-l 



X\ + i\ k ( %2 + i x 



X\ — I 



x 2 - % 



e^e^du, (2.49) 



'[0,27r)x[0,27r) 

for all k, I G Z. Let ju be another solution of the moment problem (2.1) and v be 
defined by 

u(A) = //(T _1 (A)), A G *B([0, 2tt) x [0, 2tt)). (2.50) 
By relation (2.49) we obtain that 



e ikv e a +dv= / e ik<p e il ^du, k,leZ. 

[0,27r)x[0,27r) </ [0,27r) X [0,27r) 



(2.5i; 



By the Weierstrass theorem we can approximate Lp m and ip n , for some fixed 
m, n G Z + , by trigonometric polynomials P k (tp) and R k (if)), respectively: 

max |</? m -P fc (^)| < -, max - RM)\ < -, keN. (2.52) 



Then 



/ <p m i(, n dv - [ P k {^)R k {i,)dv 

J[0,27r)x[0,27r) J [0,2tt) x [0,2tt) 



(^ m - p k {<p))i> n dv + / p^)(f - HjtCV'))^ 

[0,27r)x[0,27r) J [0,2tt) x [0,2tt) 

< max |^ m |T^([0,27r)) + max |P fc (y>)||i/([0, 2tt)) 

V>S[0,27r) A; ¥?€[0,2tt) «; 



< max |^ m |-^([0,27r)) 

^G[0,2tt) ' k 



1 \ 1 

- + max |w m | -i/( 0, 2tt)) -»■ 0, 

k ¥3£[0,2tt) Ik 



as fc — > oo. In the same manner we get 



tp^du- / P k (<p)R k W)du 

[0,27r)x[0,27r) i[0,27r) X [0,2tt) 



as k — > oo. Hence, we conclude that 



->0, 



ip m ^/j n du= <p m ij n dv, m,n(EZ + . (2.53) 

[0,27r)x[0,27r) J[0,2w) X [0,2tt) 

Since the two-dimensional moment problem on a rectangular has a unique solu- 
tion, we get v — v and \i — Ji. □ 



3. AN ALGORITHM TOWARDS SOLVING THE TWO-DIMENSIONAL MOMENT 

PROBLEM. 

As a first application of our results on the extended two-dimensional moment 
problem we get the following theorem. 
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Theorem 3.1. Let the two-dimensional moment problem (1.1) be given. The 
moment problem has a solution if and only if there exists a sequence of complex 
numbers {u m ^,i- >n ,r,t\{m,k,i;n,r,t)&i) which satisfies conditions (2.2), (2.8)-(2.11) and 

W m ,o,0;n,0,0 — s m,nj 771, TtGZ + . (3-1) 

The proof is obvious and left to the reader. 

Let the two-dimensional moment problem (1.1) be given. As it is well known 
(and can be checked in the same manner as for the relation (2.2)) the necessary 
condition for its solvability is the following: 

^ ^ C^m,n0^m' ,n' S-m+m' ,n+n' ^ 0, (3-2) 

for arbitrary complex coefficients a m)n , where all but finite number of ct m ^ n are 
zeros. 

We assume that the condition (3.2) holds. Repeating arguments of the proof of 
Theorem 2.1 we can state that there exist a Hilbert space "Ho and a sequence 

{h m ,n}m,n&+ SUcll that 

(h m , n , h m i , n ')u Q — s m+m> ,n+n> , TU, U, m! , n' G Z + . (3.3) 

Consider the following Hilbert space: 

H = Ho®(@u)j, (3.4) 

where Hj are arbitrary one-dimensional Hilbert spaces, j G N. We shall call it 
the model space for the two-dimensional moment problem. 

Introduce an arbitrary indexation in the set Q' by the unique positive integer 
index j: 

j G N !->■ w = w(j) = (m, k, I; n, r, t)(j) G Q' . (3.5) 

Suppose that the two-dimensional moment problem (1.1) has a solution [i. 
Consider the space L 2 and the following subsets in L 2 : 



L^fi — Lin{a;™a;2}m,nez + , H^,o — L^q. (3-6) 

We denote 

y m ,n ■= %T X 2> m,neZ + . (3.7) 

Notice that 

(?/m,n) 2/m',n')i? ^ m+m' ,n+n' i (3-8) 

for all m,n,m' ,n' G Z + . We shall also use the notations from (2. 35), (2. 36). 
Define the following numbers 

u m ,k,i;n,r,t ■= / x™(x 1 + i) k (x 1 -i) l X2(x 2 + i) r (x2-i) t dfj l , 



(m,k,l;n,r,t)E£l. (3.9) 
For these numbers conditions (2.2) hold and repeating arguments after the rela- 
tion (2.42) we construct a Hilbert space H and a sequence of elements {x m ,k,i;n,r,t}(m,k,i;n,r,t)en, 
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in H, such that relation (2.3) holds. We introduce the operator W as after (2.43). 
The operator W maps onto H. Set 

Hq = span{x mi o,o;n,o,o}m,nez + Q H. (3.10) 

Let us construct a sequence of Hilbert spaces Hj , j G N, in the following way. 
Step 1. We set 

fi = x w (i) ~ Ph x w (i), Hi = span{/i} 3 (3.11) 

where w(-) is the indexation in the set Q'. 
Step r, with r > 2. We set 

f r — x w ( r ) — -Pff e(®i<t<r-i-fft) a ' l(, ( r )' — span{/ r }. (3-12) 
Then we get a representation 

H = H ®(^H^j. (3.13) 

Observe that Hj is either a one-dimensional Hilbert space or Hj = {0}. We 
denote 

A M = G N : ^ ^ {0}}, A; = N\A„. (3.14) 

Then 

H = H ®\Q)H j \. (3.15) 
\ieA„ / 

We shall construct a unitary operator U which maps H onto the following sub- 
space of the model space %: 

H = H ® I 0^ I <ZH. (3.16) 

Choose an arbitrary element 

x = ^,71^^1,0,0^,0,0 + ^1\71\ — ' (3-17) 

m,n£Z + jeA M 1 1 /j II 

with a m)n , /3j G C. Set 

£/x = ^2 a m , n h rri)n +'^2 fyej, (3.18) 
where ej G "Hj, ||ej||^ = 1, are chosen arbitrarily. 

Let us check that this definition is correct. Suppose that x has another represen- 
tation: 

x= ^2 a m,n x m,0,0',n,0,0 + /^j]j"r4j — ) (3.19) 
m,neZ + jeA M ll-'jH^ 
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with 5 m) „, /3j G C. By orthogonality we have (3j = /3j, j G N. Then 

= {<Xm,n ~ a m,n) x m,0,0;n,0,0 

m,n£2,+ 

^ (^m,n C^m,n) (S^m 1 ,n' ^m',n') {hm,ni ^ , m',n')'H 



m,n,m'n 



n 



Thus, the operator [/ is defined correctly. lixEH and 



/i 



11/ 



where a m , n ,/3j G C, then 

(x,x)jj = O min a m / i ,i'(s mi o i 0;n,0,0) :(; m',0,0;n',0,o)if + 

^ ^ <-^m,n'-^m' ,n'{h"m,n-> hm' ,n , ')'H ^ ^ PjPj 

= (Ux,Ux) H . 

By continuity we extend 17 to a unitary operator which maps onto Then 
the operator £/W is a unitary operator which maps onto "H. We could define 
this operator directly, but we prefer to underline an abstract structure of the 
corresponding spaces and this maybe explains where the model space comes from. 
We set 

h m , k ,i;n,r,t ■= UWy mA i. nir>t , (m, k, I; n, r, t) G ft. (3.20) 

Observe that 

h m ,o,o;n,o,o = h m , n , m,n G Z + ; UWH^q = H . (3.21) 

Since 



i J S A M : j<r 



(3.22) 



then 



^»(r) £%e| Hj \ , r G N. (3.23) 

Observe that {y m ,k,i;n,r,t}(m,k,i;n,r,t)€n satisfy relations (2.16)-(2.19) (with y instead 
of x). Therefore {h mt k,i- n ,r,t} (m,k,i;n,r,t)en satisfy relations, as well. Notice that 

(h"m,k,l;n,r,t, h m > fi' ,1' ;ri ,r' ,t') H = (ym,k,l;n,r,ti Um' ,k' ,V ;ri ,r' ,t') L 2 
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= u m+m' ,k+V ,l+k';n+ri ,r+# ,t+r' j (3.24) 

for all (m, k, I; n, r, t) , (m 1 , k',l';n',r', t') £ f2. 
Choose an arbitrary r £ N. By (3.12) we may write 

X w { r ) = fr + ^H ffi(ei< 1 <r-iHt) Xll '( r ) 

= IIMI|i71+ E Af^ + ^o, (3.25) 

tGA M : l<t<r-l H^ll 

where ft £ C, w e # - By (3.20) and (3.18) we get 

h w{r) = Ux w{r) = \\f r \\e r + ^2 Pt e t + wo, (3-26) 

teA M : l<t<r-l 

where wq = Uuq £ "Ho- Therefore 

(/i w ( r ), e r ) > 0; (3.27) 

and 

K{r) e Ho © Ut\<* {h w(r) ,e r ) = & f r = & r £ A^. (3.28) 

yeA M : i<t<r-i / 

In particular, we may write 

r £ A M 4» (^(r), e r ) > 0, r £ N. (3.29) 

Theorem 3.2. Lei t/je two-dimensional moment problem (1.1) be given and con- 
dition (3.2) holds. Choose an arbitrary model space H with a sequence {/i m , n }m,nez + ; 
satisfying (3.3) and fix it. The moment problem has a solution if and only if there 
exists a sequence {h mt k,i- n ,r,t}(m,k,i;n,r,t)£n, in H such that the following conditions 
hold: 

1) hm,0,0;n,0,0 = ^m,n; 711,71 E Z +) ' 

2) h w (r) £ "H © (0j £ A: o<j<r^)» and ( h w(r),e r ) > 0, r £ N, for some 
subset AcN. 

3) The sequence {h m ^,i- n ,r,t\ { m ,k,i,n,r,t)£_n satisfies conditions (2.16)-(2.19) (with 
h instead of x). 

4) There exists a complex function <p(m, k, I; n, r, t), (m, k, I] n, r, t) £ Q, such 
that 

{h m ,k,l;n,r,U h m > <k > ,1' ;n> ,r> ,t>) H = f( m + m' , k + I', I + fc'j 71 + 7l' , T + t' , t + r') , 

for all (m, k, I; n, r, t), (m', k', l'\ n', r', t') £ Q. 

Proof. The necessity of conditions l)-4) for the solvability of the two-dimensional 
moment problem was established before the statement of the Theorem. 
Let conditions 1),3),4) be satisfied. Consider the extended two-dimensional mo- 
ment problem (2.1) with 

u m ,k,l;n,r,t ■= ip(m, k,l;n,r,t), (m, k, I] n, r, t) £ f2, (3.30) 
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where <p is from the condition 4). Then 

^ ^ a rn,k,l;n,r,tOim' ,k',l';n' ,r' ,t' u m+m' ,k+l' ,l+k';n+n' ,r+t> ,t+r' 

(m,k,l;n,r,t),(m' ,k' ,V ;ri ,r' ,t')efi 

= ^ ] <^m,k,l;n,r,t a m' ,k' ,V ;n' ,r' ,t' {h m ^^-n,r,ti h m ' ,k' ,1' ;n' ,r' ,t' )u 

(m,k,l;n,r,t),{m',k',l';n',r',t')en 

2 

(m,k,l;n,r,t)sQ 

for arbitrary complex coefficients a m ^,i;n,r,t, where all but finite number of ot m ,k,kn,r,t 
are zeros. 

By conditions 3) and 4) we conclude that conditions (2.8)-(2.11) hold. By The- 
orem 2.2 we obtain that there exists a non-negative Borel measure \i in M 2 such 
that (2.1) holds. In particular, using conditions 4),1) we get 



>0, 

n 



x^x^dfi = w TO! o,o;n,o,o = ^(m, 0, 0; n, 0, 0) 

{h m ,0,0;n,0,Q, ^0,0,0;0,0,o) H = (h m ,n, hofijH = s m,n, TTl,n G Z + . 



□ 



Observe that condition 2) can be removed from the statement of Theorem 3.2. 
However, it will be used later. 

Denote a set of sequences {/i m ,ifc,Z;n,r,t}(m,A;,Z;n,r,t)er2j m H satisfying conditions 
l)-4) by X = X(l-L). As we have seen in the proof of Theorem 3.2, for an arbi- 
trary {/i m ,fc,i ;n ,r,t}(m,/!,(;n,r,i)e(i £ the unique solution of the extended two- 
dimensional moment problem with moments (3.30) gives a solution of the two- 
dimensional moment problem. Observe that all solutions of the two-dimensional 
moment problem can be constructed in this manner. Indeed, let \x be an arbitrary 
solution of the two-dimensional moment problem. Repeating arguments from 
relation (3.6) till the statement of Theorem 3.2 we may write 



xf{Xx + i) {Xl - i) X% (X 2 + i) T (x 2 -iydfl= {y m ,k,l-n,r,h yo,0,0;0,0,o)L2 

= (UWy m> k,i;n,r,t, UWy Q 

,0,0;0, 0,0)% (^m,fc,/;n,r,£; 

= tp(m, k, I; n, r, t) = u m ^ n ^ 

for all (m, k, I; n, r, t) G fl. Here the operators U,W, the sequence {h mt k t i ]n ,r,t}(m,k,i;n,r,t)en, 
the function tp(m, k, I; n, r, t) and the moments u mt k,i;n,r,t, of course, depend on the 
choice of fi. Notice that {h m> k,i ;n>r ,t} (m,k,l;n,r,t)£Q 

G X(H). 

For such constructed parameters, the measure fi is a solution of the extended 
two-dimensional moment problem considered in the proof of Theorem 3.2. Since 
the solution of this moment problem is unique, \i will be reconstructed in the 
above described manner. 

Notice that condition 4) of Theorem 3.2 is equivalent to the following condi- 
tions: 

(hm,k,l;n,r,t, hm' ,k' ,l';ri ,r' ,t')u = (h> m ,k,l;n,r,t, h m ' ,k',l';n' ,r' ,t'JH; (3.31) 
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if m, m', m, fh', n, n' G Z + , k, I, r, t, k', I', r', t' G Z: m + m' = m + m'; 

{h m ,k,l;n,r,U hm',k',l';n',r',t')H = { h m ^,l;n,r,V h m' ,k' J' ;n' ,r' ,t')n, (3.32) 

if m,m',n,n' G Z+, k,k,l,r,t,k',r,l',r',t' G Z: k + V = k +T; 

{h m ,k,l;n,r,U h m > <k ' ,1' ;n' ,r> ,t>) H = { h m ,kJ;n,r,V ^m' ,%' ,l';n' ' ,r> ,t')^ (3.33) 

if m,m',n,n' G Z + , k,l,l,r,t,k' ,k' ,V ,r' ,t' G Z: I + k' = T+k f ; 

(h mt k,l- n ,r,t, h m '±',l';n',r',t')u = (h m ,k,l;n,r,t, h m > ±' ,1' ■%< ,r' ,t')ui (3.34) 

if m, m', n, n', n, n' G Z + , k, I, r, t, k', /', r',t' G Z: n + n' = n + n'; 

(h mt k t l-n,r,t, h m r t k',l'-n',r',t')H = {h m ,k,l;n,r,ti ^m' ,k' ,V ;ri ,r' t^Hi (3.35) 

if m, m', n, n' G Z+, A;, /, r, r, £, /c', r', t',t' G Z: r + £' = r + t'; 

{h m ,k,l;n,r,U h m ' ,k> ,1' ;n' ,r> ,t') H = (^m.fc^n.JvT' h m ' ±' ,1' - n ' ,r' ,t') H i (3.36) 

if m, m', n, n' G Z + , k, I, r, t, t, k', I', r', r',t' G Z: t + r' = t + r'. 

As we can see, the solving of the two-dimensional moment problem reduces to 
a construction of the set X{T-L). Let us describe an algorithm for a construction 
of sequences from X(H). 

Let {g n }^Li be an arbitrary orthonormal basis in Hq obtained by the Gram- 
Schmidt orthogonalization procedure from the sequence {h m ,n\m,n&L + indexed by 
a unique index. 

Choose an arbitrary j G N. Let w(j) = (m,k,l;n,r,t)(j) G Q'. If we had 
constructed {^ m ,n}m,nez + G X(H), then the two-dimensional moment problem 
has a solution /i and 

d j ■= \\ h w(j)\\n = \\xT(xi + i) k ( x i ~ i) l x%{x2 + i) r (x 2 - ifWh 

= / xl m (xl + l)\x\ + l) l xl n (x 2 2 + l) r (x 2 2 + lfd/I. 

Therefore dj are bounded by some constants Mj = Mj(S) depending on the 
prescribed moments S := {s mtn } m ^z + - (Notice that e.g. (x\ + 1)' < 1, for / < 0, 
and for non-negative m, k, I; n, r, t the values of dj are determined uniquely). 
Step 0. We set 

h m fl,0;n,0,0 — h mj n, 171, 11 G Z + . (3.37) 

We check that conditions (2.12)-(2.15) (with h instead of x) and (3.31)-(3.36) are 
satisfied for /v.o,o;n,o.cb fn,n G Z . If they are not satisfied, the two-dimensional 
moment problem has no solution and we stop the algorithm. 
Step 1. We seek for h w (\) in the following form: 

oo 

^u(i) = ai >n9n + Pi ; iei, (3.38) 

n=l 

with some complex coefficients ai ;n ,/3i;i. 

Conditions (2.12)-(2.15) (with h instead of x) and (3.31)-(3.36) which include 
h w (i) and the already constructed h mt k,i- n ,r,t are equivalent to a set L\ of linear 
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equations with respect to «i ;n , n G N, and d\ = Notice that they 

depend on only by d\. Denote the set of solutions of these equations by 

S\ = {(ai- n , n G N\ d\) : equations from L\ are satisfied} . (3.39) 

Set 

Si = j(ai ;n ,n e JV;di) G S 1 : ^|a 1;n | 2 < d l5 di < A/i) . (3.40) 
Finally, we set 



n=l 



Gi = < J^ai-n^n + di - ^ \ a hn\ 2 ) ■ { a l;ni n E N; d\) G i?i > . (3.41) 



n=l \ n=l 



The case G*i = is not excluded. 

Step r, with r > 2. We seek for h w r r \ in the following form: 

oo r 

^(r) = O^nflh + A-y^, (3.42) 

n=l j=l 

with some complex coefficients a r - n ,f3 r .j. 

Conditions (2.12)-(2.15) (with h instead of x) and (3.31)-(3.36) which include 
h w M and the already constructed h mi k,i- >ni r,t are equivalent to a set L r of linear 
equations with respect to a r;n , n G N, f3 r -j, 1 < j < r — 1, and d r = ||/i«,(r)||%, 
and depending on parameters (h w m, h w (2), /Wr-i)) G G>-i . 
Notice that these linear equations depend on /3 r;J - only by d r . Denote the set of 
solutions of these equations by 

S r = {(a r]n ,n G N; f3 r .j, 1 < j < r - l;d r ; h w{1) , h w{2) , /^(r-i)) : 

(h w (l), h w (fy, ^(r-l)) G Gr-1, 

and equations from L r with parameters (h w (i), h w ^), ^iu(r-i))> are satisfied} . 

(3.43) 

Set 

S r = {(a r;n ,n G N;f3 r -j,l <j<r- 1; d r ; h w (2), /i w ( r -i)) G <5 r : 

OO 1 1 ^ 

^ K;nT + ^ \Pr-j\ 2 < d r , d r < M r I . (3.44) 
n=l j=l J 

Finally, we set 

h"w(r—l) i 



1 

oo r— 1 / oo r— 1 \ 2 

^ ^ Qr;n9n ^ ^ Pr\j&j I ^ ^ |Cr;n| ^ ^ I /^"ijl 

n=l j'=l \ n=l j'=l 



G iV; /3 r; j, 1 < j < r - 1; d r ; h w ^),h w{2 ), ^(r-i)) G SV j . (3.45) 

(The case G> = is not excluded.) 

Final step. Consider a space H of sequences 

h= {ht,h 2 ,h 3 ,...), h r eH, re N, (3.46) 
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with the norm given by 

||h|| H = sup-=||/i r || w < oo. (3.47) 

r£N V M r 

For arbitrary (hi, ...,h r ) G G r , we put into correspondence elements h G H of the 
following form 

h= (hi,...,h r ,g r+ i,g r+2 , ...) : gj G H, \\gj\\ n < s/Mj, j > r. (3.48) 

Thus, the set G r is mapped onto a set G r C H. If G r = 0, we set G r = 0. 
Observe that all elements of G r has the norm less or equal to 1. Set 

oo 

G = pi G r . (3.49) 

r=l 

If G 7^ 0, then to each (gi,g 2 ,...) G G, we put into correspondence a sequence 
^ = {h m ,k,i;n,r,t}(m,k,i;n,r,t)en such that (3.37) holds and 

h w (r) := g r , r G N. (3.50) 

We state that fj G X(H). In fact, conditions (2.12)-(2.15) (with h instead of x) 
and (3.31)-(3.36) are satisfied for h m0)0 . n)0t0 , m,n G Z +) by Step 0. If one of these 
equations include h w M with r > 1, then we choose the maximal appearing index 
r. Since (h w ^, h w ( r )) G G r , then this equation is satisfied. Condition 2) is 
satisfied by the construction. 

Thus, if G 7^ 0, then using Sj we can construct a solution of the two-dimensional 
moment problem in the described above manner. 

Theorem 3.3. Let the two-dimensional moment problem (1.1) be given and con- 
dition (3.2) holds. Choose an arbitrary model space H with a sequence {/i m ,„} mi nez + , 
satisfying (3.3) and fix it. The moment problem has a solution if and only if 
conditions (2.12)-(2.15) (with h instead of x) and (3.31)-(3.36) are satisfied for 
h m ,o,0;n,o,o '■= h m ^ n , m,n G Z +; and 

G^0, (3.51) 

where G is constructed by (3.49) according to the algorithm. 

If the latter conditions are satisfied then to each (gi,g2, •••) G G ; we put into 
correspondence a sequence = {h mt k,i; n ,r,t}(m,k,i;n,r,t)£n such that (3.37) holds and 

h w (r) ■= 9r, r G N. (3.52) 

This sequence belongs to X(l-i) and the unique solution of the extended two- 
dimensional moment problem with moments (3.30) gives a solution /i of the two- 
dimensional moment problem. Moreover, all solutions of the two-dimensional 
moment problem can be obtained in this way. 

Proof. The sufficiency of the conditions in the statement of the Theorem for 
the solvability of the two-dimensional moment problem was shown before the 
statement of the Theorem. Let us show that these conditions are necessary. 
Let /i be a solution of the two-dimensional moment problem. By Theorem 3.2 the 
set X(H) is not empty. Choose an arbitrary = {/i m ,fc,i ; n,r,f} G X(U). 

By conditions 3), 4) we see that conditions (2.12)-(2.15) (with h instead of x) and 
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(3.31)-(3.36) are satisfied. In particular, they are satisfied for /i m ,o,o;n,o,o = h mt7l , 
m,n G Z+. 

Comparing condition 2) with Steps 1 and r for r > 2, we see that 

h w(r) ) G G r , r G N. (3.53) 

Therefore elements 

(h w (i ) ,...,h w ( r) ,g r+ i,g r+2 ,...)eG r , r G N, (3.54) 
where <?j G H : Hg^H^ < y/Mj, for j > r are arbitrary. Thus, the element 

h := (h w (i),h w{2 ),K {3) , ...) G G r , r G N. (3.55) 

and 

h G p| G r = G. (3.56) 

r-eN 

Therefore G ^ 0. 

If the conditions of the Theorem are satisfied then to each g = (<7i, <?2j •••) £ 
G, we put into correspondence a sequence = f)(g) = {/j m ,fc,i in ,r,i}( m ,fc,iin,r,t)eJi 
such that (3.37) and (3.52) hold. Then Sj G X(H), as it was shown before the 
statement of the Theorem. The sequence Sj generates a solution of the extended 
two-dimensional moment problem and of the two-dimensional moment problem, 
see considerations after the proof of Theorem 3.2. 

It remains to show that all solutions of the two-dimensional moment problem can 
be obtained in this way. Since elements of X(H) generate all solutions of the two- 
dimensional moment problem (see considerations after the proof of Theorem 3.2), 
it remains to prove that 

{55(g): gEG} = X(n). (3.57) 
Denote the set on the left-hand side by X\. It was shown that X\ C X(9i). 

On the other hand, choose an arbitrary = {/i m ,/M;n,r,t} G X(H). 

Repeating the construction at the beginning of this proof we obtain that 

h := (h w (i),h w (2),h w (3), ...) G G. (3.58) 

Observe that 

£(h) (3.59) 

Therefore X(H) C X l and relation (3.57) holds. □ 

Remark 3.4. The truncated two-dimensional moment problem can be considered 
in a similar manner. Moreover, the set of indices of the known elements h m ^ n will 
be finite in this case and therefore equations in the r-th step of the algorithm 
will form finite systems of linear equations. Thus, the r-th step could be easily 
performed using computer. 

Remark 3.5. Consider the following system of r linear equations: 

/ fi \ 




(3.60) 



V ti ) 
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where A 1 = (aij)i<i<r;jeN is a given complex numerical matrix, 1 < i < r, are 
given complex numbers, and Xj, j G N, are unknown complex numbers; r G N. 

The Gauss algorithm allows to solve this system explicitly. Let us briefly 
describe this. 

Step 1. If A 1 = then the algorithm stops. Conditions of solvability and the 
set of solutions are obvious in this case. 

If A 1 ^ 0, let mi-th column of A 1 be the first non-zero column of A 1 . Inter- 
changing equations we set the non-zero element of this column in the first row 
and divide this equation by this element. Then we exclude x mi from the rest of 
equations. We get the following system: 

%m-L + a l,m 1 +l X m.\+l + a l,m 1 +2 X m 1 +2 + ••• = f\ ; (3.61) 



^mi+1 
A 2 f X m 1+ 2 



( fi\ 

fi 



(3.62) 



V/ r 2 / 

where A 2 is a given complex numerical matrix with r — 1 rows, f 2 , 1 < i < r, 
and afj, j > m\ + 1, are given complex numbers. 

Then we repeat the same construction for the linear system (3.62). After a 
finite number of steps the algorithm stops. Then we exclude x mt ,x mt ^, ■■■,x mi 
from the previous equations (t < r). 

Thus, the numbers xf j ^ m 1 ,m 2 , ■■■,mt can be chosen arbitrary such that the 
corresponding series in (3.60) converge, and x mi , x m2 , ...,x mt are defined uniquely. 
If t < r, we additionally have the solvability conditions which follow from (3.62) 
in the last step. 

Observe that if we have an infinite number of equations in (3.60), we can 
choose an increasing number of equations and then construct the intersection of 
the solution sets. 

A modified algorithm. Notice that in Theorem 3.3 the correspondence between 
the parameters set G and solutions of the two-dimensional moment problem is not 
necessarily bijective. The algorithm may be modified to make this correspondence 
one-to-one. The following modified algorithm is more complicated. If we only 
need to check the solvability or the bijection is not necessary for our purposes, 
we can use the original algorithm. 

First, condition 2) of Theorem 3.2 may be replaced by the following more precise 
condition: 

2) Set A := {r G N : (h w{r) ,e r ) > 0}. Then 

h w{r) G U © f HA, and (h w{r) ,e r ) > 0, r G N. (3.63) 

VjGA: l<j<r J 

The necessity of this condition was shown before Theorem 3.2, while condition 2) 
was not used in the proof of the sufficiency of Theorem 3.2. 
As before, we denote a set of sequences {h mt k t i ]n ,r,t}(m,k,i;n,r,t)en, in H satisfying 
conditions l)-4) of Theorem 3.2 by X = X(%). Observe that the modified X 
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is a subset of the original one. The same arguments show that the new X{H) 
generates all solutions of the two-dimensional moment problem, as well. 
Step and Step 1 of the algorithm will be the same as before. 
We set 

H k :=H ® ( Hj J , U\ := {h G H k : (h,e k ) > 0}, k G N. (3.64) 

\l<j<k / 

In the r-th step we shall proceed in the following way (r > 2). 
Choose an arbitrary (/Wi), h w (2), h w r r -i)) G G>_i. Observe that by the con- 
struction in the (r — l)-th step we have 

h w{j) G U k ~ l or h w(j) eV. k + , 1 < k < r - 1. (3.65) 

Set 

:= {s — (si, S2, s r -i) : Sj — 1 or Sj — 0, 1 < j < r — 1}; (3.66) 

and 

Hp := {(h u h 2 , /i r _!) : hj G "H^ 1 if Sj = 0; hj G ?4 if Sj = 1; 1 < j < r-1}, 

s6 5 r . (3.67) 
Observe that S r is a finite set of 2 r_1 binary numbers. By (3.65) we obtain that 

(h w (i), h w (2), /i«,( r -i)) G GV-i nHJ, for some s G S"\ (3.68) 

Set 

r r _ M := G r -i n Hp, s G ,S r . (3.69) 

Notice that 

rv-i^j n r r _i )?2 = 0, si, s 2 e s r , (3.70) 

and 

G r _i = |J r r _ 1>y . (3.71) 

Choose an arbitrary s*G S ,r . We seek for /i„>(r) in the following form: 

00 

h>w(r) ^ ^ OLr;nQn ^ ^ fir;j&j Pr\r&ri (3.72) 

n=l l<j<r— 1: Sj—1 

with some complex coefficients a r;n , P r .j\ /3 r;r > 0. 

Conditions (2.12)-(2.15) (with h instead of x) and (3.31)-(3.36) which include 
h w t r ) and the already constructed fo m & j ;n r ,t are equivalent to a set L r (s) of linear 
equations with respect to a r;n , n G N, f3 r] j, 1 < j < r — 1 : Sj = 1, and 
(i r = || h W (y) ||%) and depending on parameters (h w (i),h w (2), h w ^ r -i)) G r r _i^ . 
Notice that these linear equations depend on /3 r - t j only by d r . Denote the set of 
solutions of these equations by 

S r (s) = {(a r - n ,n G N; P r -j, 1 < j < r - 1 : Sj = 1; d r ; h w ( X ),h w (2), ^(V-i)) : 

(/iu,(l), /lio(2)> •••) ^io(r-l)) £ r r _i ; ^, 

and equations from L r (s) with parameters (^(i), h w (2), ^(r— i)), are satisfied} . 

(3.73) 
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Set 

S r (s) = {(a r . n , n G N; (3 r - :j , 1 < j < r - 1 : Sj = 1; d r ; h w{1)j h w ( 2 ), /i w ( r -i)) G SV 

oo 

5^ l«r;n| 2 + ^ |A-yf < <*r, d r < M r I . (3.74) 

n=l l<j<r-l: Sj=l 

Finally, we set 

G r (s) = { h w (2), ^w(r-l)) 



OO 



n=l l<j<r-l: s,=l \ n=l l<j'<r-l: Sj=l 



r;j I 



12 ) e r 



(a r . n ,n G N; f3 r;j , l< j <r - I: Sj = l;d r ; h w ^,h w ^ 2 ), ■-, ^(r-i)) G <SV(s)} . 

(3.75) 

The case G r (s) = is not excluded. 
We set 

G r := (J G r (s). (3.76) 

seS r 

The final step is the same as for the original algorithm. Thus, we obtain the set 
G. 

To each g = (gi,g2, •••) G G, we put into correspondence a sequence .fj = fj(g) = 
{h m ,k,i;n,r,t}{m,k,i;n,r,t)eu such that (3. 37), (3. 50) hold. We state that i}(g) G X(U). 
Observe that the set G r in (3.76) is a subset of G r for the original algorithm. 
Therefore the set G is a subset of G for the original algorithm. Thus, .fj(g) 
belongs to the old X{%). To show that it belongs to the modified X(H) it 
remains to verify (3.63). Observe that (h w m, h w (2), h w r r )) G G r (s), for some 
s G S r . The condition (h w r r ), e r ) > follows from the construction of h w t r ) in the 
r-th step. Set 

A ( r \ — I I 1 < J < r - 1 : Sj- = 1} U r, if (h w{r) ,e r ) > , , 

y W-| {l<j<r-l: S , = l}, if(^ (r) ,e r ) = • 

By (3.72) we get 

^ (r) G Ho © I ^ j . (3.78) 

yeAo(r) / 

Thus, it remains to verify that A (r) = {j G A : 1 < j : < r} =: A(r). But for 
1 < j < r — 1, conditions = 1 and (h w (j\, ej) > are equivalent. Consequently, 
we obtain f)(g) G X(H). 

Theorem 3.3 remains true if we replace words "according to the algorithm" by 
the words "according to the modified algorithm", and add the following sentence: 
" The correspondence between elements of G and solutions of the two-dimensional 
moment problem is bijective". Let us check this last assertion (and the rest of 
the proof is similar). 
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The correspondence between G and X(l-L) is obviously bijective. Let Sjj = 
{ h L,k,i;n,r,t}{m,k,i;n,r,t)an G X (H) , j = 1,2, be different: ^ fi 2 - They pro- 
duce solutions yt,i and f/,% of the two-dimensional moment problem, respectively. 
Suppose that [i\ = /i 2 = /i. Recall that fij is constructed as a solution of the corre- 
sponding extended two-dimensional moment problem with moments u 3 mkl . nrt = 
(Pj(m, k, l;n,r,t), j = 1,2 (see the proof of Theorem 3.2). Here ipj(m, k, I; n, r, t) is 
from Condition 4) for S)j, j = 1, 2. Therefore <pi(m, k, I; n, r, t) = ^(m, k, I; n, r, t). 
By condition 4) of Theorem 3.2 this means that 

{^ l m,k,l;n,r,ti^ l m',k',l';n',r',t')'H = (^m,k,l;n,r,t: ^m',k',l'-,n',r',t')'Hj (3.79) 

for all (m, k, I; n, r, t) , (m', k', l';n',r', t') G f2. 
Choose the minimal r, r G N, such that 

h l w{r) ± h 2 w{r) . (3.80) 

By (3.37), (3. 79) we obtain 

Pn Kv(r) — P-Ho^wir) =: ^0- (3.81) 

By condition (3.63) we may write 

h i(r) = h + ^2 7rtf e j + lr;r e r, ll-j & C, 7* P > 0; 

l<j<r-l: (hl uy e 3 )>0 

h 2 w{r) = h + Yl Ir^j + ^r, ieC, 7r 2 ;r >0. (3.82) 

l<j<r-l: (hl u) ,e 3 )>0 

Since h 1 ^ = h 2 w{j) =: h w(j) , 1 < j < r - 1, we get 

{j: 1 < j < r - 1, (^ (J) , e,) > 0} = {j : 1 < j < r - 1, y) , e .) > 0} =: A. 

(3.83) 

Therefore 

K(r) = h o + Yl ftj e i + ^V- %3 G C ' ^ir > 0, a = 1, 2. (3.84) 

Suppose that there exists j G A such that 7I7 7^ 7 r 2 ;J -. Let j be the minimal such 
index j. Since j £ A, we get 

Go := (h w (j )i e j ) = (Kt>(j )i e jo) > 0' a = 1,2, 

and 

^C?o) = Oo e io + "jo-ii «j -i G (3.85) 

Then 

e jo = — (h w{jo) - % _i); (3.86) 

and 

Tryo = (Ku(r)i e jo) = ^=(^C(r)> hw(jo) — u jo-l) 
Vo 

1 1 
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= =(h a w[r) ,h w{jo) )-=(h + 7rVi' u io-0- (3-87) 

By (3.79) and our assumption about jo we obtain 7* - o = 7 2 j - This contradiction 
means that 7^ = 7 2 .-, Vj G A. Therefore 

^C(r) = h + 7r;r e r, ^ £ W~\ % r > 0, a = 1, 2. (3.88) 

Observe that 

ll^ W l| 2 =II^H 2 +l7, a ;r | 2 ' a=1 ' 2 - 
By (3.79) we conclude that jj. r = 7 2 r . Therefore = We obtained 

a contradiction with (3.80). The proof of the last assertion for the modified 
Theorem 3.3 is complete. 

4. ON A CONNECTION WITH THE COMPLEX MOMENT PROBLEM. 

In this section we shall analyze the complex moment problem: to find a non- 
negative Borel measure a in the complex plane such that 

z m z n da = a min , m,n G Z + , (4.1) 

where {a min } mi „ e z + is a prescribed sequence of complex numbers. 
Recall the canonical identification of C with IR 2 : 

z = x\ + X21, X\ = Re z, X2 = Imz, z G C, (x%, x 2 ) G IR 2 . (4.2) 

Let a be a solution of the complex moment problem (4.1). The measure a, viewed 
as a measure in M 2 , we shall denote by \i a . Then 

= — — — V V CFC?(-l) n ~ j / z k+] z m ~ k+n ^da 

K I k=0 j=0 JV 
^ m n 

= 2 m (2i) n EX^ -1 )™ J CrC" a fc+j,m-fc+n-j, (4.3) 
where C£ = Then 



c 



z m z n da — (xt + ix 2 ) m (xi - ix 2 ) n d(i (7 
Jr 2 

m n „ 

J^J^ ^ ^- 1 ^' 1 / ^l + '(^2) m " r+n "^/ia 
„ — n i— n */M 2 



r=0 /=0 

m n 



\ \ 1 nm/m/ 1 \n— I •m— r+n— i „ 

/ J / J W °Z l — X J 4 *r-M,m-r+n-Zi 



r=0 Z=0 
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and therefore 

m n 

a m ,n = VV C™Cn-l) n - l * m - r+n - l s r+ i,m-r + n-i, m,neZ + , (4.4) 



r=0 1=0 



where 



^(-l) n - j CJ?C?a k+jt m -k + n-j, m,neZ + . (4.5) 



2 m (2i 

v i k=0 j=0 

Since \i a is a solution of the two-dimensional moment problem, then conditions 
of Theorem 3.3 hold. 

Theorem 4.1. Let the complex moment problem (4-1) be given. This problem 
has a solution if an only if conditions of Theorem 3.3 and (4-4) with s m ^ n defined 
by (4-5) hold. 

Proof. It remains to prove the sufficiency. Suppose that for the complex moment 
problem (4.1) conditions of Theorem 3.3 and (4.4) with s mn defined by (4.5) hold. 
By Theorem 3.3 we obtain that there exists a solution \i of the two-dimensional 
moment problem with moments s m ,n- 

The measure /i, viewed as a measure in C, we shall denote by cr M . Then 



z m z n da fl = (xi + ix 2 ) m (x 1 - ix 2 ) n d\i 
c Jr 2 

m n „ 

= ^^C r m C7(-l) n "' / x[ +l (ix 2 ) m - r+n - l d/i 

r=0 1=0 ^ R2 

m n 

\ \ ^ /~im(~m( i\n~l-m—r+n-l„ 

— / / 1 h r+l,m-r+n-l ~ a m,m 

r=0 1=0 

where the last equality follows from (4.4). □ 

Theorem 4.2. Let the complex moment problem (4-1) be given and conditions 
of Theorem 3.3 and (4-4) with s m ^ n defined by (4-5) hold. Let ^ be a set of all 
solutions of the complex moment problem (4-1) and $ be a set of all solutions of 
the two-dimensional moment problem (1.1) with s mj „ defined by (4-5). Then 

^ = {a M : /j£<£}. (4.6) 

Therefore all solutions of the moment problem (4-1) are described by Theorem 3.3. 

Proof. The proof is straightforward. □ 

Of course, this Theorem holds for the modified version of Theorem 3.3, as well. 
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